Effective interactions in the stripe state of the two-dimensional Hubbard model 
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We investigate the structure of the pairing potential in the stripe phase of the two-dimensional 
Hubbard model. Based on the random phase approximation we discuss in detail the interactions 
in the charge- and spin channel and compare our calculations with related considerations in com- 
mensurate antiferromagnets. Our main finding is that due to the incommensurate charge-density 
wave formation the exchange of collective modes in the charge channel is significantly enhanced 
with respect to the spin bag approach whereas due to the inhomogeneous charge distribution the 
coupling to transverse spin fluctuations tends to be suppressed. 
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I. INTRODUCTION 



The existence of electronic inhomogeneities is now a 
well established experLmpiiital fact in underdoped high 
Tc materials (see e.g.ljl3). Especially in lanthanum 
cuprates these inhomogeneities can develop in the form 
of antiphase domain walls where quasi one-dimensional 
charge stripes are separated by hole-free antiferromag- 
netically (AF) ordered regions. It is worth to note that 
stripe textures have been predicted as stableJHartree- 
Fock (HF) saddle-points of the Hubbard modelDEl before 
they were found experimentally in cuprates and nicke- 
lates. Stripe textures in high-Tc materials have been 
first detected in Nd-doped lanthanum cuprate by neu- 
tron scattering experiments where due to the occurence 
of a low temperature tetragonal phase (LTT) both in- 
commensurate AF and charge order are pinnedQ. The 
relevance of structural distortions for the formation of 
antiphase domain walUs has t"fpi recently studied in 
Ref.t3. In Nd-free compoundsMO these measurements 
have up to now only detected the magnetic part of the 
scattering which remarkably displays the same incom- 
mensurability than in the Nd-doped system. In addition 
evidence for incommensurate charge fluctuations comes 
from the analysis of optical phonon measurements- in 
YBa2Cu306.6t3 and neutron diffraction experimentsEil in 
La2-a;Sra;Cu04 (LSCO). Taking the stripes as an exper- 
imental fact the question arises wether they have some- 
thing to do with high-Tc superconductivity or if they are 
just some 'strange byproduct' of the strong correlations 
which are at work in these materials. In fact there are 
many experiments which support the point of view that 
stripe correlations could give rise to high transition tem- 
peratures. First , it has been shown that the stripe spac- 
ing in LSCO compounds is inversely proportional to the 
transition temperature up to optimal dopingtll. A close 
connection between transition temperature and stripe in- 
commensurability has also been found in YBCOliS where 
interestingly there seems to appear a discontinuous jump 
of incommensurability from 1/6 at optimal doping to 1/8 
for the Tc=60K underdoped compound. Second, the 
occurence of domain walls can be iHterpreted in terms 
of a quantum critical point (QCP)li3 located near opti- 



mal doping Xopt which is a quite appealing concept in 
order to account for the anomalous normal state trans- 
port of high-Tc materials around this particular doping. 
As a consequence the QCP scenario allows for a natural 
subdivision of the phase diagram into ovcrdoped metal- 
lic and-i;he underdoped pseudogap phase. It has been 
arguedO that the QCP can be deduced from the ampli- 
tude of the incommensurate magnetic peaks vansishing 
near Xopt . This coincides with experiments where super- 
conducting order has been suppressed by strong pulsed 
magnetic fieldsES and which have revealed an underlying 
metal insulator transition at about the same concentra- 
tion in agreement with the QCP concept. Temperature 
dependent measuxements of the charge and spin^sector in 
Nd-doped LSCOcI as well as NQR experimentsEl further 
support the idea that it is the charge rather than the 
spins which is reponsible for the stripe state instability. 

The possibility of d-wave pairing near a incommepj 
surate CDW instability has been investigated in Ref.Ej 
within a BCS-type scheme. This analysis was restricted 
to the overdoped and optimally doped region of the 
phase diagram where according to the QCP scenario 
strong ICDW fluctuations are present but without sym- 
metry breaking of the translational invariance. In the 
underdoped system the opening of a ICDW gap natu- 
rally leads to a suppression of the superconducting order 
parameter due to a reduction ia. the density of states 
around the Fermi level (see e.gH3,El). However, in ad- 
dition the ICDW fluctuations are strongly modified in 
the symmetry-broken state but to our knowledge no de- 
tailed analysis of the corresponding pairing interaction 
has been performed yet. The present paper is dedicated 
to explore the effective electron-electron interactions in 
the underdoped regime of the phase diagram thus start- 
ing with a ground state which exhibits long-range incom- 
mensurate stripe order. Our considerations are based on 
the Hartree-Fock (HF) decoupled two-dimensional Hub- 
bard model where we obtain the domain wall structure 
via a self-consistent iteration of the on-site charge and 
spin expectation values. From the formal point of view 
our investigations can lie viewed as a generalization of 
the spin-bag approachEj to incommensurate antiferro- 
magnetism. The basic idea of the spin-bag mechansim is 
that a hole doped into the commensurate antiferromag- 
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net locally perturbes the AF order parameter and con- 
sequently two holes gain energy when they share a com- 
mon deformation i^ag) in analogy to a standard lattice 
bipolaron. In RefEj this scenario has been worked out by 
considering the pairing in the longitudinal spin channel of 
a commensurate antiferromagnet which indeed turns out 
to be attractive for small momentum transfers whereas it 
becomes repulsive for large q ^ (tt, tt). The charge chan- 
nel plays a minor role for pairing in the spin-bag approach 
since it does not exhibit an instability and renormaliM*^ 
tion effects are rather weak. It has been argued in Ref.E3 
that the repulsion in the transverse spin channel at large 
q may be suppressed by the coherence factors entering 
the vertex of the pairing interaction. However, as pointed 
out by Frenkel and Hankec3 the vanishing of the vertex 
is cancelled by the Goldstone pole in the spin susceptibil- 
ity thus leading to a finite repulsive contribution in the 
transverse channel for momenta q ~ (vr, tt) comparable 
to that for the amplitude fluctuations. 

We want to emphasize that our HF approach to the 
stripe phase fails to account for some of the experimen- 
tal facts mentioned above but on the other hand allows 
for a simple description of some basic features regard- 
ing the electron-electron interaction in a gwstem posess- 
ing stripe order. First, it is well knownEJ that a HF 
decoupling of the one-band Hubbard model cannot prop- 
erly describe the observed stripe charge structure in the 
cuprates which is characterized by one hole for every sec- 
ond unit cell along the domain wall ('half-filled'). In con- 
trast HF calculations result in completely filled stripes 
thus predicting twice the observed stripe periodicity and 
it has turned out that correlations beyond HF theory 
have to be includcctjO in order to account for the mea- 
sured stripe structure within the Hubbard model. 

Thus in our investigations we will restrict to excita- 
tions perpendicular to the domain walls which should 
be less affected by the actual stripe filling than excita- 
tions along the walls which of course crucially depend on 
the doping of the stripes. Second, it has been shown in 
Ref.E3 that the HF stripe instability is due to an insta- 
bility in the spin channel and the corresponding charge 
modulation is generated as a higher harmonic of the lon- 
gitudinal spin scattering. Thus in our model the spins 
are the driving force behind the domain wall formation 
in contrast to what is observed in Nd-doped LSCO. An 
alternative approach where antiphase domain walls are 
enslaved by a frustrated phase seperation scenario has 
been discussed incj. However, since we restrict on zero 
temperature we believe that within the limitations dis- 
cussed above our approach can capture some features of 
the electron-electron interaction in the stripe state of un- 
derdoped cuprates. We note that a similar study for the 
calculation of correlation functions has been performed in 
Ref.E3 where the authors focussed on the renormalization 
of the spin velocity in the stripe state. 

The paper is organized as follows. In Sec. II we out- 
line the formalism of our approach starting from HF the- 
ory and the subsequent derivation of effective hamiltoni- 



ans for the residual RPA-type interactions in the stripe 
phase. In Sec. Ill we present results for these interac- 
tions focussing on the different contributions from ver- 
tex functions and susceptibilities respectively. We finally 
conclude our discussion in Sec. IV. 



II. FORMALISM 
A. Mean-fleld theory 

We start with the Hubbard hamiltonian 

= ^ efcC^Cfc,,, + ^ ^ (1) 

where c^^]^ destroys (creates) an electron in the state k 

and energy efe, pq^a = J2k A+q.a'^k.a denotes the density 
operator. 

For one elecron per site the ground state shows two- 
sublattice commensurate antiferroniagnetic order. It has 
been shown by SchuliZI using the RPA-type Stoner crite- 
rion that for a small number of holes the system posesses 
a magnetic instability at an incommensurate modulation 
wavevector (7r,7r) -I- qg. In principle the corresponding 
(mean-field) symmetry-broken state could be realized ei- 
ther by a modulation of the transverse spin degrees of 
freedom i.e. a spiral phaseEj or the corresponding mod- 
ulation of jthe longitudinal magnetization. From general 
arguments and also from unrestricted Hartree-Fock cal- 
culations (see e.g.cil) one finds that usually the symmetry 
breaking occurs in the longitudinal channel of the mag- 
netic system. As a consequence this leads to a strong 
coupling between charge- and spin degrees of freedom re- 
sulting in the formation of antiphase domain walls where 
the charge is distributed along the phase boundaries. It 
should be noted that the inclusion of the vacuum renor- 
malization of the effective interaction may stabilize the 
commensurate antiferromagnet at small dopiiigsLJ. 

Within unrestricted HF theory (see e.g.Ej) one finds 
that for small values of the on-site interaction U the 
stripes are oriented along the verticals and whereas for 
large U they run along the (l,l)-direction. In the fol- 
lowing we consider a rather small value of U close to the 
instability in order to realize the vertical stripe phase. 

Moreover we restrict ourselves to an energy dispersion 
arising from nearest neighbor hopping = —2t[cos + 
cos ky) with bandwidth B = 8t ^ AeV . For doping 
5 — 0.13 (measured from half-filling) one finds from 
the Stoner criterion the occurence of the magnetic in- 
stability at Q* = 7r(l ± 1/8, 1) and Q* = 7r(l, 1 ± 1/8) 
when the Hubbard interaction exceeds the critical value 
Uc ~ 0.9ey. In principle these four equivalent wave vec- 
tors allow either for an unidirectional stripe structure or 
for a two-dimensional grid pattern of antiphase domain 
walls. In the spin glass phase of LSCO where the stripes 
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run along the diagonals neutron scattering provides coa- 
vincing evidence for a one-dimensional modulationEZI, 
however, in the superconducting regime this paint re- 
mains unclear although from general arguments^ unidi- 
rectional scattering seems to be favored. 

Restricting ourselves to the (longitudinal) domain wall 
structure the Hubbard interaction can be decoupled via 

P-q-aPq,a ~ ^ {{P-q)Pq + {Pq)P-q - {P~q)Pq " {Pq) P-q) 



A(rr) 



1 {{Pq){P-q) - {Pq){P-q)) 



(2) 



where {p^} = Y.aiPq,<y) and {pq) = Y.a'^(Pq.'y) denote 
the charge and spin densities respectively. Although the 
Stoner criterion signals the instability at a single wave 
vector Q* — (tt, it) +qs the anharmonicity of the domain 
wall structure in the symmetry broken phase naturaU* 
incorporates higher harmonics also. It is well knownEfl 
that the topological nature of the stripe phase implies 
the relation qc — 2qs between the basic charge (qc) and 
spin (qs) modulations. 

In order to illustrate the role of these harmonics in the 
charge- and spin channel respectively, we have sketched 
in Fig. 1 the charge (xi^)) ^nd spin order {A{x)) pa- 
rameter for a simple rectangular stripe pattern in the 
periodicity interval [0, /]. One can easily convince oneself 
that the corresponding Fourier coefficients Xk and Afc are 
given by 

Xk = -^cos(fc7r/2)cos(— ^) 

Afc = — ^ sin(fc7r/2) sin(— ^) (3) 

KIT I 

where ^ denotes the width of the magnetic domains. 
Thus the magnetic scattering only involves odd harmon- 
ics whereas the charge scattering is due to the even har- 
monics. This fact holds in general when both magnetic 
and charge stripes are symmetric under reflection with 
respect to their central axis. As a result one can describe 
the modulation of charge and spin densities in terms of 
the single wave vector Q* 

(p(i-)>=2 (/OpQt)cos(pQ*r), 

p{even)—2 

{p{r))=2 (MpQOsin(pQ*r) (4) 

p(odd) — l 

where the commensurability M of the stripe is de- 
fined as the number of scattering events in the same 
direction needed to return to the equivalent state, i.e. 
|k-|-MQ*) — |k). Note that within this definition a large 
commensurability M not necessarily implies a scattering 
vector Q* very distant from (tt, tt). 



FIG.l. Illustrative plot for the cliarge (xi^)) a-nd spin order 
(A(a;)) parameter of a rectangular stripe pattern. 

For the vertical stripe orientations along the x- 
direction considered here we can write Q* = + 
Q{M — 2),1) where the step function is defined with 
0(0) = 0. This guarantees the correct behavior for 
commensurate AF (M=2) i.e. Q* = (7r,7r). Note that 
within these definitions the commensurability M corre- 
sponds to the period of the stripe structure (in units of 
the lattice constant) as is indicated in Fig. 2a. Fur- 
thermore we are now able to construct the reduced Bril- 
louin zone (BZ) which is spanned by vectors Q* and 
Q-L = i2+MeiM^2)P+AfA -M,2+MeiM-2)). In order 
to illustrate the new zone scheme we plot in Fig. 2b the 
reduced BZ for the case M = 16 which will be studied be- 
low. Note that the reduced BZ is not only due to charge 
order (the corresponding vertically striped BZ is also in- 
dicated in Fig. 2b) but is rotated and divided in half due 
to the presence of (modulated) antifcrromagnctism. 

Introducing reduced zone notation for the electron op- 
erators 

C„(fc, cr) = Ck^ + {n-l)Qi,ky + (n-l)Ql,a (5) 

we thus obtain the following hamiltonian 
H=Y <^n{k)c\{k,a)cn{k,a) 

n,k,a 

+ XI [X2pci+2p{k,<j)cn{k,a) + h.c.\ 
+ X! cr['^pci+p{k,cr)cn{k,a) + h.c] 

p,7i,k,a 

+ ^5](|A,p-|X2/) (6) 

p 

where we have introduced the charge and spin or- 
der parameters as X2p — {U/2){p2pQt) and Ap — 
-(C//2)(/.pQ0. ^ 

The hamiltonian eq. (|6|) can be diagonalized via 

M 

Cn{k, cr) = ^ Anm{k, (j)fm{k, a) (7) 
m— 1 
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FIG. 2 a) Sketch of a vertical stripe pattern where the shaded 
areas correspond to the charged domain waUs. The commen- 
surability M times the lattice constant a defines the period- 
icity of the structure, b) The large square represents the full 
Brillouin zone whereas the dark shaded rectangle depicts the 
reduced Brillouin zone for incommensurability M = 16. The 
latter can be constructed from the reduced zone for charge or- 
der only (light shaded rectangle) which is rotated and divided 
in half due to antiferromagnetism. Since also the AF order 
is modulated the boundaries of the dark shaded rectangle are 
not parallel to those of the reduced zone for commensurate 
antiferromagnetism (dashed line). 

and from the symmetry of H it turns out that the 
transformation matrices obey the relations Anm{k) = 
Anm{k,(T) = (-I)"v4„„(fc,-CT) and A„„(fc,o-) = 

We finally obtain the hamiltonian in diagonal repre- 
sentation as 

k.ij n 



where En{k) denotes the dispersion of the n-th sub- 
band in the reduced BZ (dark shaded rectangle in Fig. 
2b). The order parameters have to be determined self- 
consistently. 

In Fig. 3 we show a cut of the bandstructure through 
the reduced zone and the density of states for a stripe 
structure with commensurability M=16 corresponding 
to the scattering vector Q* as defined above. As can 
be seen the spectrum is composed of 16 quasi one- 
dimensional subbands. Since the choosen carrier con- 
centration 5 = 0.13 slightly exceeds the integer fraction 
5 = 1/8 the chemical potential is located within but near 
the top edge of the 7th subband. 

The following sections deal with the RPA fluctuations 
around the HF stripe state and it is quite instructive to 
consider the functional dependence of the ground state 
energy from the spin- and charge order parameters which 
is shown in Fig. 4. In fact, since RPA excitations can be 
derived from the expansion around the HF saddle-point 
in the density fluctuations (see e.g.Lj) one can already 
qualitatively deduce from Fig. 4 the possible instabili- 
ties of the system. Note that for the E{Ai) curve the 
charge order parameter is adjusted selfconsistently while 
for the E{x2) curve the spin order parameter is taken at 
the saddle point. The {{H),Ai) curve displays the stan- 
dard behavior of a second-order phase transition towards 
the SDW ground state with minimum at Ai = O.lAeV. 
On the other hand it turns out that the saddle-point of 
((i?), X2) corresponds to a maximum in the energy-order 
parameter space which is also obvious from the decou- 
pling eq. (||)c3. This simple picture already demonstrates 
that in the present framework charge order is completely 
due to the SDW instability. 




(0,0) ^ f(|,l) 



0.2 0.4 0.6 0.8 

D0S(1/eV) 



;(iA,p 



\X2p\ 



(8) 



FIG. 3 Band structure and density of states for a stripe array 
with commensurability M=16 and U/t=1.3. The wave vec- 
tor cut of the bandstructure is from (0, 0) along Q* (see Fig. 
2b). The dashed line indicates the position of the chemical 
potential. 
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FIG. 4 Free energy as a function of the SDW order parameter 
Ai. The inset displays the free energy as a function of the 
CDW order parameter X2- 

Moreover, note that for X2 — > the {{H),X2) curve 
has a finite derivative which signals the absence of crit- 
ical fluctuations in the charge channel for the homoge- 
neous system in contrast to models where the instability 
mechanism is due to frustrated phase separationca'ta. 



B. Susceptibilities for general commensurability 

We start from the following definitions for charge and 
spin-density correlations functions 

Xm'^m' = -^{Tpq+{.m-l)Q\a-{t)P-q'-{m'-l)Q'',cr'{Q)), 
Amm 



(9) 



where the wave vectors q are within the reduced zone. 
Since the mean-field hamiltonian eq. (H) commutes with 
Sf the excitations in the longitudinal spin and charge 
channel are decoupled from the transverse channel where 
the solution of the RPA approximated Dyson equation 
for the latter case leads to 

Xt^-^M,^) = [1 - uZ\:r{q.Lo)]-\'^\:-{q.co) (10) 
The bare transverse stripe susceptibilities read as 



k s t 



nt{k)[l -ns{k + q)] 



ns{k + q)[l~ nt{k)] 



uj + Es{k + q) - Et{k) - irj u; + Es{k + q) - Et{k) + irj 
where we have introduced the vertex functions: 



^ttmik - 



^q, fc) = E A*j^Jk + g)A„,i(fc)(-l)"-l 



rr,„(fc,fc + g) = (r+„(fc + q,fc))* (12) 

and Ant{k) is defined in Eq. (7|). Here the index 
J\f = J\f{n,'m,k,q) defines the number of the BZ which 
contains the unreduced state \k + q) and k + q denotes 
the corresponding reduced wave vector. 

We now turn to the evaluation of the correlation func- 
tion in the longitudinal channel. In this case we start 
from the RPA decoupled Dyson equation for Xm'rji' ^) 

x™m'(9,^) = Xw''(g,^) + t^Xw'^lg, w)x„?',';^(<7,c^) 

(13) 

where 

'"(9''^) ^ ^J2^stnAkl^q,k)T^sm'ik,kl^q) 

k st 

nt{k)[l - ns{k + q)] 



Am rtv 



ns{k + q)[l - nt{k)] 



uj + Es{k + q) - Et{k) - ir] w + Es{k + q) - Et{k) + irj 
and the vertex functions are given by 

^1tm{k + q. k)^Y. A*^Jk~q)A^,t{k)a"+^ 



^{k,k + q) = {T-^^{k + q,k)r 



(15) 



Since we have to project Xm^n'i^^^) ^^^^ ^^^^ charge 
and longitudinal spin sector we introduce the new matri- 
ces 



1 



pp 



pz 

mm Amm' 



A- 



1 



Amm' Amm' 
Amm' Amm' 



where 



A. 



(16) 



(17) 



Within these definitions we can rewrite eq. (^_3|) as a 
2M X 2M matrix equation for the charge and spin density 
correlations respectively 



s = [1 - [/s°r^]-i5° 

and the matrices T, and S read as 



(18) 



/ Kli Ki2 • • • \ 

K21 K22 



V • 



/r^ 



V • 



(19) 



Q Obviously both charge and longitudinal spin excita- 
tioris are strongly coupled except in the half-filled system 
where the ground state corresponds to a commensurate 
antiferromagnet . 



(14) 
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C. Effective interactions for general 
commensurability 

Having determined the collective excitations in the 
striped phase we are now in the position to study the in- 
teractions of two holes within the same subband with op- 
posite spin and momenta in the spin- and charge channel 
respectively. Note that the susceptibilities defined in eq. 

acquire non-diagonal matrix elements in the stripe 
state, which implies the possibility of creating Cooper 
pairs with a finite momentum (which is a multiple of Q*). 
However, for simplicity we will consider only scattering of 
zero momentum Cooper pairs and therefore neglect non- 
diagonal contributions of the susceptibilities (m ^ m') to 
the interaction. Additionally we will cast these interac- 
tions in the form of effective hamiltonians thus restricting 
on the static limit with an appropriate frequency cutoff. 

We find in the orientational spin-fluctuation channel 



rL„(fc,g)ri,„(fc,g) = (-i)('"-i)5, 



(25) 



Now remember that the index m in Xmm(q) labels 
the multiple of Q* which has to be added to the re- 
duced transfered momentum q, i.e. the transfered mo- 
mentum in the full BZ is q = q -I- (m — 1)Q*. Thus 
an odd (even) index m corresponds to enhanced charge 
(spin) susceptibilities which therefore in general medi- 
ate attractive (repulsive) interactions respectively. In 
addition upon assuming that the charge (spin) suscep- 
tibility has some finite contribution iar odd (even) m 
only, one recovers the result of Ref£3 that the effec- 
tive interaction in these channels becomes proportional 
to the commensurability of the charge modulation, i.e. 
E™rL„(/c,9)r^,„(/c,g)=±Af/2. 



N ^ 



^{k + q, k) r+ „(-A: - g, -A:)x™m('7) 



III. RESULTS 



X ft{k + q,-cj)ft{-k-q,a)ft{k,a)ft{-k,-a) (20) 
in the charge-fluctuation channel 

^'^ = E ^'stra{k + q. k) r,T™(-fc^ 9, ~k)x'„'.^{q) 

m st 

X /+(fcTq,-a)/+(-fc^g,a)/t(fc,a)/t(-fc,-a) (21) 



and in the amplitude spin-fluctuation channel 

stmik + q,k) T-^,^{-k - q, ~k)x"nrni<l) 



X ftik + q,-<j)f+i-k-q,a)Mk,a)M~k,-a). (22) 

In order to illustrate the behavior of the vertex func- 
tions in the various channels let us consider the case of 
well developed stripe order (i.e. large charge- and spin 
order parameters). In this limit we can as a first ap- 
proximation neglect the kinetic energy in eq. (^ and the 
transformation cq. (^ takes the form 



Ant{k) 



1 r 27r , 

^=exph — nt\. 
/M m 



(23) 



As a result the vertex contribution in the transverse 
spin channel is given by 



r;,,„(fc,g)r+„(fc,g) = i~i)^"'-'^s,,^M 



(24) 



and thus the coupling to spin-flip scattering vanishes for 
intraband transitions. 

On the other hand the intraband coupling in the charge 
and longitudinal spin channel survivesCj and can be ex- 
pressed in the form 



We now apply the formalism developed above to the 
calculation of the effective interaction within the HF 
stripe state of the two-dimensional Hubbard model. We 
restrict ourselves to the domain wall considered in Sec. 
Ila, i.e. commensurability M = 16 for which the total 
scattering vector is Q* = {tt^tt + Fig. 5 displays the 
different RPA susceptibilities for ui — and uj = 25meV 
respectively. Note that the q-scans for the magnetic sus- 
ceptibilities are along the {qx,qy — tt) axis whereas the 
charge susceptibility is shown along {q^ ,qy — 0) . The 
structure in is dominated by the Goldstone pole at 
Q' due to the breaking of spin-rotational invariance in the 
stripe state. As can be seen from Fig. 6 the mode splits 
with increasing frequency and the two braiiches result 
from the defolding of the reduced BZ (Ref.t3). Moreover 
beyond the gap energy (roughly corresponding to Ai) 
the magnetic excitations rapidly lose intensity. 

The longitudinal spin and charge channel do not dis- 
play a physical pole aX uj — Q. However, a mode like fea- 
ture is present at Q* for the amplitude fluctuations and 
at 2(5* for the charge susceptibility which reflects the en- 
hanced susceptibility of the stripe state to perturbations 
with the respective wave numbers. Since both charge and 
spin are coupled within the RPA formalism the respec- 
tive excitations follow the same dispersion in (g, ijj)-space 
(Fig. 6) but since the charge 'mode' appears as a second 
harmonic to the longitudinal spin excitation its intensity 
is significantly lower. It is interesting to observe that 
at the energy (« Ai) where the longitudinal spin fluc- 
tuations become commensurate the charge fluctuations 
show the periodicity of the static spin modulation (i.e. 
qs) which again is a consequence of the defolding of the 
reduced BZ. Note that this energy does not correspond to 
the maximum in the intensity, which appears for slightly 
lower frequencies. 
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(0,7r) q (7r,7r 

50 I ■ ■ 

"> Lo — Q meV 



(0, tt) q (7r,7r) 
Lo — 2h meV 



(0,^) 



(tt, tt) 



w = me\ 




(0,7r) q (7r,7r 
;', w = 25 meV 




(0,0) q (^,0) (0,0) g (^,0) 



FIG. 5 The momentum dependence of the RPA susceptibili- 
ties XrpA'XrpajX'rpa ™ ^he stripe state for two frequences 
oj = (left column) and lo = 25meV (right column); Solid: 
real part. Dashed: imaginary part; first and second lines rep- 
resent the scans along path (0, tt) — (tt, vr) , third line represents 
the scan along path (0, 0) — (tt, 0). 

Finally, in order to determine the efRcieny of the pair- 
ing in the various channels we have to calculate the ver- 
tex functions which determine the coupling of the charge 
carriers to the corresponding correlation functions. Gen- 
erally speaking the vertex functions determine the over- 
all sign of the effective interaction in the Brillouin zone. 
Let us briefly illustrate this effect in the limit of a com- 
mensurate spin-density wave (M = 2) where the ver- 
tex functions can be calculated analytically. It is then 
shown below that the qualitative behavior of the cou- 
pling still holds for higher commensurability. Note that 
the pairing interaction mediated by amplitude fluctua- 
tions of AF order constitutes the basis of the so called 
"spin bag mechanism to high-Tc superconductivity" as 
proposed by Schrieffer et al. in Ref.BS. In this case the 
antiferromagnetic scattering Q = (tt, tt) doubles the unit 
cell and we denote conduction and valence band by the 
superscripts c and v respectively. The corresponding en- 
ergies are given by E'^^'" = ±Ek — ±-y/A2"+^ and A is 
the SDW order parameter. 

Restricting ourselves on the 'hole doped' case the rele- 
vant coupling functions can be written within our nota- 
tion as 



Obviously there is a difference in sign between small 
q- (m = 1) and large q- (m = 2) scattering. Since in an 
antiferromagnet the magnetic susceptibilities are peaked 
at (tt, tt) the attractive interaction for small momentum 
transfer (as anticipated in the spin-bag mechanism) is 
accoiapanied by a large repulsion at q ~ Q (see e.g. 
Ref.H,Ea for a more detailed analysis of this point). 

Let us now consider the case of commensurability 
M=16. As already mentioned the considered hole doping 
6 = 0.13 is slightly above the 'magic' concentration 1/8 
so that the Fermi level is located within the 7th subband 
in the reduced BZ. 

In Fig. 7 we have plotted the corresponding vertex 
functions for intraband scattering which for convenience 
are defolded along (qxtO) and (qxtTt). It turns out that 
also in this case the coupling around odd multiples of Q* 
is negativ whereas it is positive around the even harmon- 
ics. Thus also in the stripe phase the large q scatter- 
ing which is dominated by the magnetic fluctuations is 
repulsive whereas attraction predominantly arises from 
the coupling to charge fluctuations. Naturally for stripes 
the latter are much more pronounced in comparison to 
the spin-bag model. With regard to the effective hamil- 
tonian eqs. ( pO[p^ ) the cutoff of the pairing potential is 
an important parameter in determining the interaction 
strength. From the right panel of Fig. 5 we conclude 
that for the relevant wave- vectors near the 'modes' the 
cutoff energy approximately is of the order of the charge 
order parameter X2- 

Besides the enhancement of charge fluctuation medi- 
ated pairing in the stripe phase the contribution of spin- 
flip scattering to the interaction is also signiflcantly differ- 
ent from the spin-bag approach. First note that the cou- 
pling to the transverse susceptibility for large momenta 
(upper right panel of Fig. (7)) is by two orders of magni- 
tude smaller than the corresponding coupling to charge 
and longitudinal spin fluctuations. In addition it vanishes 
at the momentum of the Goldstone pole i.e. q — Q*. For 
the commensurate M — 2 case this can be seen directly 

(tt, tt) (note that in our 
in eq. (|27 



from eq. (27) in the limit q 

reduced-potation this is equivalent to q 

In RefB this fact led to some confusion concerning the 
relevance of the transverse fluctuations within the-|Spin- 
bag approach. However, it has been shown in Ref.u that 
the vanishing of the vertex is cancelled by the gapless 
Goldstone pole of (Q) and as a result the interaction 
in the spin-flip channel around (tt, tt) becomes compara- 
ble to the pairing from the amplitude fluctuations in the 
spin-bag scenario. 




FIG. 6 The momentum and frequency dependences of the 
imaginary parts of the RPA susceptibihties Xrpa'X'rpa^x'rpa 
in the stripe state; right column shows the intensity con- 
tour along (0,7r) — (vr,7r) for Im{x^pj^),I'ni(xjiPA) ^-nd along 
(0,0)-(7r,0) for /m(x?fp^).. 

In order to study the coupling to transverse fluctua- 
tions for commensurability M = 16 we have performed 
a numerical expansion of the vertex and the transverse 
susceptibility around q = Q*. Contrary to the low com- 
mensurability case one finds that the contribution from 
the transverse channel to the pairing is only « 5% of that 
of the spin amplitude fluctuations. We consider this as 
a consequence of the spatial separation of charge carriers 
and spin fluctuations in the stripe model. 



IV. CONCLUSION 

To summarize, we have calculated the effective inter- 
actions arising from the RPA fluctuations around the 
HF stripe phase of the two-dimensional Hubbard model. 
Since our approach is a generalization of Ref.E3 to higher 
commensurabilities the general features resemble closely 
those of the spin-bag mechanism. However, significant 
differences arise due to the inhomogeneous charge dis- 
tribution in the stripe state where the holes are con- 
fined to the domain walls separating antiphase AF re- 
gions. As a result the charge correlations are enhanced 
for wave-vectors corresponding to the stripe periodicity 
resulting in a much more pronounced small q attrac- 
tive pairing as compared to the spin-bag approach where 



the system is translationally invariant with respect to 
the charge degrees of freedom. Due to the fact that 
HF theory generates the charge order as a second har- 
monic of the longitudinal spin modulation the pairing is 
still dominated by the large q repulsion due to magnetic 
fluctuations. However, in models where the stripe in- 
stability is driven by the charge an analogous approach 
should yield a comparable interaction strength in spin- 
and charge channels respectively. Thus depending on the 
repulsive contribution due to magnetic scattering either 
a d-wave or anisotropic s-wave SC order parameter min- 
imizes the free energy of the system. Obviously non of 
these order parameters will be really realized in a purely 
one-dimensional stripe phase as considered in the present 
paper, however, our investigations can be considered as 
complementary to Ref.c^ where the coexistence between 
incommensurate antiferromagnetism and anisotropic su- 
perconductivity has been investigated. In this context 
the effective interaction derived above may provide a mi- 
croscopic origin_for the superconducting correlations con- 
sidered in Ref.Ej. Anyhow, in order to have true d-wave 
order one has to anticipate a scenario where the stripes 
fluctuate between the x- and y- directions in order to 
restore the two-dimensional symmetry. In principle one 
could symmetrize the pairing interaction derived above 
and use it in a BCS-type approach in order to evaluate 
the superconducting gap. Alternatively, one could con- 
sider the interactions arising from the fluctuations of an 
eggbox-type domain wall structure as in Ref.a in order 
to mimic orientationaL stripe fluctuations. Furthermore 
it has been proposedcj that domains of stripes running 
along the crystalographic directions can couple via the 
twin boundaries which also leads to 'symmetric' d-wave 
superconductivity. 




(0,0) q (^,0) (0,^) q (^,^) 




FIG. 7 The vertex functions l'7rrni^^l)^77mi^^l) ^^'^ 
r77„(fc, g)r^7'^(fe, g) for two points in the reduced BZ:fc = 
(— 7r/1.6, — 7r/2.3), solid curve; k = (0,0), dashed curve; The 
dotted line indicates the position of Q*. 
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However, within the obvious Hmitations of our ap- 
proach which we have aheady discussed in the introduc- 
tion a quantitative evaluation of the superconducting gap 
would not be very instructive and beyond the scope of 
the present paper. Here our purpose is to simply outline 
the general features of stripe fluctuation induced pairing 
and to eventually supplement considerations within more 
sophisticated approaches. 
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